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Abstract — Recent experimental and theoretical (both improved semiclassical and fully 
quantum-mechanical) line broadening calculations for Bill and NeVII An = transi- 
tions with n = 2 and 3, respectively, are discussed. The disagreements by about a factor 
of 2 between the fully quantum-mechanically calculated and both measured and semi- 
classically calculated widths can be explained in terms of violations of validity criteria 
for the semiclassical calculations and nonthermal Doppler effects. Only the quantum 

calculations allow a clear separation of elastic and inelastic scattering contributions 
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ion scattering experiments. Additional independent determinations of line widths for 
multiply-ionized, nonhydrogenic ions are called for, but meanwhile caution should be 
exercised in the use of corresponding semiclassically calculated widths, e.g., in opacity 
calculations or for plasma density diagnostics. 



1. INTRODUCTION 



Most spectral lines observed in laboratory experiments or lines emanating from stellar 
^S^ , atmospheres are, as far as their Stark broadening is concerned, in the category of isolated 



lines. Such lines, per definition, have widths, not to mention shifts, which are smaller than 
separations between the perturbed and the relevant perturbing levels. Relevant in this 



context are especially levels connected to upper or lower levels of the transition by electric 
dipole matrix elements, i.e., levels which would contribute to the Stark shift of the line 
in an external electric field. Furthermore, if these widths are small enough, such that the 
required timescale for the Fourier transform of the line profile (in angular frequency units) 
is large compared with the duration of important electron-ion collisions, then the impact 
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approximation to the general theory can be used to calculate widths (and shifts) of the 
resulting Lorentzian profiles. 

Most important is generally the (full) width, 7, which according to Baranger 1 can be 
written as 



1 = N e vF{v) 
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(0,v)~f e (0,v)\ 2 dn\dv, (1) 



in terms of cross sections a uu i and aw for clcctron-collisional inelastic transition rates from 
upper (u) and lower (I.) levels into the perturbing levels u' and £' ', respectively. These cross 
sections are functions of the initial electron velocity v, and the aw-product is averaged using 
the electron velocity distribution function F(v) as weight function. The ensuing inelastic 
contribution to the line broadening, often called lifetime broadening, is as expected, e.g., 
from the original work of Lorentz. 2 However, the last (elastic) scattering term is less intuitive. 
It involves the difference of elastic scattering amplitudes / on upper and lower levels. These 
quantities depend on velocity v and scattering angle 0, and were it not for the cross term 
~fufe ~ fufe ' ^ ne integrals over v and O would simply give the sum of total elastic scattering 
rates on upper and lower levels. However, this cross term can be very important, especially 
for An = transitions, and cause strong cancellations. This could already have been 
inferred from the semiclassical, adiabatic phase shift limit, 3 ' which had been applied also 
for the case of monopolc-dipolc interactions 6, (quadratic Stark effect) at the time of the first 
realistic Stark broadening calculations of isolated lines from neutral atoms 6 and ions 8,9 in 
low charge states. These calculations all used the semiclassical approximation and were only 
done to the lowest nonvanishing order in the interaction, the phase-shift calculations being 
a notable exception to this. 6 ' 7,10 Results of such calculations, and a dispersion relation 11 
between widths and shifts, were actually used to optimize the choice of minimum impact 
parameters and strong-collision terms in perturbative semiclassical calculations of widths 
and shifts. 12 For neutral and singly-ionized atoms this procedure is surprisingly accurate, 
especially for the widths. 13 In other words, remaining errors are much smaller than the 
strong-collision contributions given in Tables IVa and V of Rcf. 12, which are rather large. 
For lines from multiply-ionized atoms, both measurements and calculations have been 
much more difficult and are still relatively rare. The experimental problems lie in the 
development of suitable, i.e., optically thin and homogeneous, plasma sources capable of 
obtaining the high densities and temperatures required, and in the application of density and 



temperature diagnostics. Most successful in this regard has been the gas-liner pinch work at 
Bochum; 14 comparisons 15 of some of these results with linear-pinch-discharge experiments 16 
are also possible and interesting. As to calculations, various versions of a semiempirical 
method 17 had been rather successful 13 in predicting widths of many lines from low charge 
states. However, along isoelectronic sequences, and for An = 0, n — 3 transitions in 
multiply-ionized Li-like, 18 Be-like 19 and B-like 20 ions, significant disagreements were found, 
measured widths being larger than most calculated widths for the higher members of the 
isoelectronic sequences. For An = 0, n = 2 transitions, measured widths for (Li-like) B III 21 
and, perhaps, also for C IV 22 disagree with most semiclassical and semiempirical estimates 
in a similar way. 

The semiempirical estimates were based on Eq. (1) in the sense that the effective Gaunt 
factor approximation 23,24 was used for the cross sections of dipole-allowed collision-induced 
transitions, e.g., from the upper level of the line, with an extrapolation to bclow-thrcshold 
energies (velocities) to account for elastic scattering from ions in the upper state of the line. 
At least in Ref. 17, the lower level width calculated in the same way was simply added, 
i.e., the interference term, — /„// — fuft > mentioned above was neglected, as were collision- 
induced transitions caused by other than dipolc (A = 1) interactions, i.e., monopole (A = 0), 
quadrupole (A = 2), etc., interactions. The major difference between these semiempirical 
and previous versions of semiclassical calculations is that in the latter the equivalents of 
the Gaunt factors are actually calculated using classical path integrals of the dipole in- 
teraction, supplemented by estimates of strong-collision contributions. Recently proposed 
improvements to this procedure will be discussed in Sec. 3 after a more detailed discussion 
of experiments in Sec. 2. Fully quantum-mechanical calculations are reviewed in Sec. 4, 
followed by a summary and some recommendations. 

2. EXPERIMENTS 

The stringent requirements for accurate measurements of Stark broadening parameters 
have been summarized before, in connection with the first critical data evaluation for isolated 
ion lines, 25 but are difficult to reconcile with the high densities and temperatures necessary 
for multiply-ionized atoms. Moreover, at that time, electron densities remained well under 
10 18 cm -3 and temperatures below 5 eV. Additional critically reviewed data 26-27 (see also 
Ref. 28 for an updated bibliography) range to higher densities, by about an order of mag- 
nitude, and to temperatures as high as 20 cV. The high-aspect ratio, cylindrical, gas-liner 
(gas-puff) z-pinch 14 was designed to facilitate measurements of multiply-ionized ions under 
optimal conditions, which can be well and independently (of Stark broadening) measured. 



In this gas-liner pinch, a driver gas, hydrogen or helium, is injected axially as a concentric 
hollow gas cylinder near the wall of the discharge vessel, while a test gas is injected in 
controlled (small) amounts near the axis. In this manner one ideally obtains optically thin 
line emission from a nearly homogeneous, albeit transient, plasma near the discharge axis. 
Great care is taken to have sufficient pre-ionization for the driver gas to facilitate uniform 
breakdown of the main discharge and to avoid, hopefully, instabilities during the implosion 
of this low-aspect ratio gas-puff pinch. Typical implosion velocities are about 10 cm/sec, 
corresponding to Reynolds numbers over 10 4 , suggesting the development of hydrodynamic 
turbulence. 29 (Magnetic field effects are not likely to inhibit such turbulence because of the 
relatively large ion-ion collision frequency, namely v w 10 9 sec -1 for ~ 100 cV proton-proton 
collisions, compared to proton cyclotron frequencies of u C i w 10 8 scc _1 at typical fields of 
IT.) The timescale for development and decay of the turbulence 30 is r « £/Av, if I is a 
characteristic length, say 1 cm, and Av « 10 6 cm/sec a typical spatial difference of velocities. 
This suggests r « lfisec, somewhat less than the time to peak compression, but longer than 
the duration of significant line emission. 

While the driver gas may therefore be highly turbulent during the time of interest, the 
heavier test-gas ions have collisional mean- free-paths of > 0.1cm against collisions with 
protons; they may therefore average over several eddies, especially during the implosion. 
In any case, the test gas ions can have significant nonthermal velocities, which could be of 
the same order as thermal velocities of the driver gas and therefore much larger than their 
thermal velocities. 

It is not clear to what extent the otherwise very powerful (collective) Thomson-scattering 
diagnostic 31,32 is capable of distinguishing between thermal and nonthermal (test gas) ion 
motions. (This would be different in case of plasma turbulence associated with nonlinear 
plasma waves.) Moreover, it usually measures temperatures and electron density only in a 
rather small fraction of the emitting volume, which may not be entirely representative of 
the average conditions along the line of sight for the spectroscopic measurements. Only very 
recently, 22 ' 33 with the help of a 2-dimensional detector array, have scattering spectra also 
been taken radially-imaged to check the homogeneity in regard to density and temperature 
and to verify that, e.g., argon as test-gas was restricted to a small, near-axial region. Axial 
imaging of emission spectra with krypton as test gas 3 as function of test-gas concentra- 
tion suggests a stable plasma column for concentrations below 1%, but inhomogeneitics at 
higher levels, indicating macroscopic instabilities. However, these observations cannot be 
interpreted as evidence against fully developed, fine-scale turbulence in case of small test-gas 



concentrations as used for line profile measurements. In principle, the effective tempera- 
ture for the Doppler width of lines of test-gas ion could be determined from the width 
of the impurity peak in the scattering spectrum, 31 ' 32 but the corresponding error is fairly 
large at the small concentrations needed to avoid self-absorption of strong lines. Finally, to 
conclude 19 from the absence of bulkshifts of scattering spectra (i.e., from the absence of net 
velocities of the emitting plasma volume) that there is no additional (to thermal Doppler 
broadening) Doppler broadening of the lines from test-gas ions, is premature. However, 
these observations do indicate that any turbulent eddies are smaller than 1 mm. 

For narrow lines, the possibility of systematic errors in Stark widths due to underesti- 
mates of the Doppler broadening 29 evidently remains. It is made even more plausible, in 
the case of CIV An = 0, n = 2 transitions, 22 by the observation of a Doppler splitting 
of the lines corresponding to less than or about 4 x 10 6 cm/sec, i.e., about half the implo- 
sion velocity. So far it is not understood why these streaming velocities are not completely 
randomized by collisions or turbulence. 

3. SEMICLASSICAL CALCULATIONS 

Stimulated by the various disagreements between measured widths and semiclassical 
calculations, which were already discussed in Sec. 1, two improvements to these semiclas- 
sical calculations have been proposed. The first of these was an attempt by Alexiou 35 
to improve the estimates for strong-collision impact parameter and strong-collision term 
by insisting that unitarity of the path integrals of the electron-ion interaction energy be 
preserved also during the interaction, not only over the entire collision as in previous semi- 
classical calculations. This distinction is important for hyperbolic perturber paths, because 
of the compensation between effects during approach and separation occurring especially for 
dipole interactions. (These cancellations had been noticed already in the context of electron- 
collisional broadening of ionized-helium lines. 36 ) Strong-collision impact parameters in the 
usual calculations are of the order 12 

h n 2 X n 2 

Pmin = 77 = o ly i \^l 

mv Z 27T Z 
A being the deBroglic wavelength of the perturbing electron and Z the (spectroscopic) 
charge, e.g., Z = 3 for Bill, etc. This first improvement is therefore probably not realistic 
unless n 2 /2irZ £ 1, which is not fulfilled for the lines showing large disagreements with 
previous semiclassical calculations. The 2tt in Eq. (2) was apparently omitted in Rcf. 35, as 
was the effect of curvature on the distance of closest approach in estimating the validity of 
the long-range dipole approximation. 



The second improvement, 37 by the same author, was to replace the lowest order, dipole 
interactions, perturbation theory by a numerical solution of the time-dependent equations 
to all orders, and to include long-range quadrupole interactions. This nonperturbative 
approach no longer requires the unitarity check, but of course still a minimum impact 
parameter, etc., to avoid errors from breakdown of the scmiclassical approximation and 
of the long-range approximation to the exact electron-ion Coulomb interactions. Both of 
these errors were again severely underestimated 15 in the same manner as discussed in the 
preceding paragraph. 

In spite of these potentially serious problems, the nonperturbative scmiclassical 
calculations 37 ' 38 give widths, e.g., for the NeVII 2s3s-2s3p singlet and triplet lines which 
agree within reasonable errors with gas-liner experiments, 38 whereas previous calculations 
yielded smaller widths by as much as a factor of 2. A similar discrepancy is found if com- 
parison is made with a fully quantum-mechanical close-coupling (CC) calculation, 39 albeit 
only for the rather similar 3s-3p transition in NeVIII. 

As will be discussed in the following sections, such large disagreements may be the rule 
rather than the exception for strong isolated lines from multiply-ionized atoms, i.e., for lines 
for which n 2 jZ is near or below 1. 

4. QUANTUM-MECHANICAL CALCULATIONS 

Fully quantum-mechanical, CC, calculations of electron impact broadening were pub- 
lished almost 30 years ago for resonance lines of Mgll 40 and Call, 41 and a few years later 
also of Be II. 42 It was already noticed at that time that the close-coupling (CC) calculations 
gave widths which were less than half the measured width for Be II, and less than a third 
of the semiclassically calculated width. Except for the impact approximation, 1 the colliding 
electron is treated in these calculations exactly and on an equal footing with the (active) 
bound electrons, with which it interacts via the complete Coulomb interaction rather than 
only its long-range multipolc expansion. In other words, except for numerical errors from 
using, e.g., an insufficiently complete system of basis functions, insufficient number of ac- 
tive bound electrons, or insufficient energy resolution, the CC calculations should be almost 
exact. They have therefore also been continued as part of the Opacity Project with results 
published for lines from some transitions in Li- and Be-like ions 39 and for An = 0, n = 2 
transitions in C III. 43 For the latter lines, elastic scattering contributions to the width are 
particularly important, providing 20-55% of the line widths. 

Meanwhile, completely quantum-mechanical calculations were also made for the electron 
broadening of H-like 44 ' 45 and He-like 46,47 multiply-charged ions (overlapping lines). Here 
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the emphasis was on line spectra from Z ~ 6 ions, encouraging the use of the distorted- 
wave 44,47 approximation, thus significantly reducing the numerical effort in comparison with 
CC-calculations. As a matter of fact, it turned out that the additional phase shifts even 
for L = partial waves remained small so that no unitarization or higher order calculation 
was required. (L is the angular momentum, in units of h, of the scattered electron.) The 
L value corresponding to the minimum impact parameter in semiclassical calculations, see 

Eq. (2), is 

n 2 

and a corresponding effective strong collision term was determined in Rcf. 44 by a fit to 
the sum over partial waves of the expression for the line width. (Remember, however, that 
n 2 /Z r^ 2ir and therefore L m [ n ~ 2ir is required for the semiclassical approximation.) All 
of these calculations used the full Coulomb interaction, and (red) shifts were calculated as 
well in Refs. 45-47. They are mostly due to the A = 0, penetrating monopolc interactions, 
which for the Is and Is 2 lower states are also important for the interference between elastic 
scattering terms. 

Returning to isolated lines, the discrepancies between measured widths and nonpertur- 
batively, semiclassically, calculated widths on one side, and the fully quantum-mechanical 
calculations just mentioned, have encouraged some new calculations 29 ' 48 using the converg- 
ing close-coupling (CCC) method. 49 ~ 51 This method is a standard close-coupling approach, 
except that discrete and continuum target states are obtained by diagonalizing the Hamil- 
tonian in a large orthogonal Lagucrre basis. The coupled equations are formulated in mo- 
mentum space; the convergence can therefore be easily tested by increasing the basis size. 
This method has been shown to give very similar results 52 for inelastic and elastic electron- 
ion cross sections to the RMPS (R-matrix with pscudo states) method, 53 but avoids the 
difficulties associated with the oscillatory behavior of wave functions in coordinate space. 

The first use of the CCC method for line broadening calculations was for the case of the 
B III 2s-2p (Li-like) resonance doublet measured 21 on the gas-liner pinch. Special attention 
was paid to the elastic scattering term in Eq. (1), which decreases much faster with electron 
energy E than the l/E decay of the 2s and 2p (non-Coulomb) elastic cross sections. The 
cancellation between upper and lower level elastic scattering is substantial; simply using 
the sum of elastic cross sections would lead to an overestimate of the elastic scattering 
contribution to the line width by a factor of about 6 at T e = 10 eV. However, because of 
this cancellation and the somewhat erratic energy dependence, there is more uncertainty 
in the elastic than in the inelastic contributions. Collision-induced 2s-2p excitation and 
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deexcitation transitions give the major contribution to the line width, and CCC, RMPS and 
Coulomb-Born with exchange 54 (CBE) calculations all give very similar cross sections, the 
latter being typically 20% larger than the two strong-coupling results. Most importantly, 
over 90% of the total cross section comes from L £ 6, i.e., from the nonclassical region. 
Finally, inelastic scattering associated with An > 1 transitions, estimated using the CBE 
approximation, was found to contribute only 5% to the calculated width, vs. typically 35 
and 60% from elastic and An = inelastic scattering, respectively. Moreover, the CCC 
width was even smaller, by less or about 10%, than the first i?-matrix method results 39 in 
the temperature range of interest (5-10 cV). 

Further evidence for the severe discrepancies mentioned above is obtained by a compar- 
ison of new CCC calculations 48 for the 2s3s-2s3p singlet and triplet lines of NcVII (Be- like) 
with experiments 19 ' 38 and nonperturbative semiclassical calculations. 37 In this case, two 
bound electrons are actively involved in the scattering process, while the Is 2 innershell 
electrons can be considered as part of a frozen core. (This was verified 48 by comparisons 
between the full Hartree-Fock and frozen-core Hartree-Fock calculations using the Cowan 
code. 55 ) In analogy to the Bill case, 3s-3p inelastic and super-elastic transition rates are 
again a major contribution to the line width, followed by 3p-3d transition rates. For these 
dipole (A = 1) cross sections about 50% of the total cross section arises from L > 9 for 
the singlet s-p and triplet p-d transitions, and from L > 7 for the triplet s-p and singlet 
p-d transitions. The L J 2-7T criterion is thus fulfilled only marginally. Also, the distance 
of closest approach for, e.g., an L = 5 classical electron is about equal to the boundstate 
radius, 15,48 which causes the long-range interaction used in the semiclassical calculations to 
fail. The semiclassical results are therefore again rather questionable. Note also that CCC 
and CBE cross sections for these large dipole cross sections are within 10% of each other, 
and that dipole cross sections for 2s3i?-2p3^ transitions contribute about 10%, 48 although 
they were ignored in semiclassical calculations. The An > 1 collisional cross sections for 
excitation or deexcitation of the 3£ electron are at the percent level, but were also included 
in Rcf. 48. 

The 2s3s-2s3d quadrupolc (A = 2) transition rates contribute only about 3% to the total 
quantum-mechanical inelastic transition rate, 48 also at T e = 20 eV, in part because of the 
relatively large threshold energies (~ 9cV). This is in contrast to the semiclassical result, 37 
to which the quadrupolc (A = 2) channel contributes 15% of the total width. Because this 
total width is about twice the quantum-mechanical result, the A = 2 cross sections actually 
differ by an order of magnitude. This is not very surprising, because in this case L i 5 



partial waves are most important; semiclassical and long-range interaction approximations 
are both clearly inappropriate. This might have been inferred by analyzing Fig. 11 of Rcf. 
56, in which nonperturbative, semiclassical contributions to the line width are shown as a 
function of impact parameter. The excluded region here corresponds to L m 3 for the cases 
of only A = 1 and for both A = 1 and 2 interactions. 

This leaves the elastic scattering contribution, which is essentially ignored in the semi- 
classical calculations. For this contribution, only L 5 2 partial waves are important, which 
are entirely in the quantum-mechanical and short-range interaction (A = 0, etc.) regimes. 
As for B III, there is again strong cancellation of upper- and lower-level scattering ampli- 
tudes, by about an order of magnitude, for the electron energies of interest. The contribution 
to the electron collisional line widths is only about 8 and 10%, respectively, for singlet and 
triplet lines. The estimated errors in the quantum calculation, 48 using a comparison of 
results obtained from Eq. (1) and from a formulation in terms of T-matrix elements, 12 > 57 
from CCC calculations, are actually ±15%, mostly related to resonances in the electron-ion 
scattering. 

5. SUMMARY AND RECOMMENDATIONS 

It should be clear to the reader that in contrast to a previous assertion, 58 there is still 
no convincing convergence between measurements and calculations of electron-collisional 
widths of strong isolated lines from multiply-ionized atoms, and for that matter, not even for 
the 2s-2p resonance line of Bell. 42 As far as the measurements are concerned, the observed 
trends along isoelectronic sequences are particularly striking. This observed slower decay 
with spectroscopic charge Z than predicted by most calculations is especially pronounced 
for lines which are relatively insensitive to electron broadening in comparison with higher 
n-lines. This slower decay occurs for n = 3 lines and probably also for n — 2 lines. It 
is, therefore, suggested that some of the other line broadening mechanisms compete with 
the broadening by electrons. Of these, Stark broadening by ions comes to mind, but has 
been shown theoretically 12 not to be sufficient to explain discrepancies by a factor of about 
2. (Note that there is no question concerning the validity of semiclassical calculations for 
perturbing ions.) The proposed additional Doppler broadening due to turbulent or other 
nonthermal motions of the emitting ions provides a more likely explanation. It had been 
invoked before to explain the unexpectedly large widths of singlet and triplet 2s-2p lines of 
CV in laser-produced plasmas. 59 

Hopefully the uncertainties in the interpretation of measurements of relatively small 
Stark widths will be removed, e.g., by measurements on other plasma light sources, by 



precision measurements of the impurity peak in collective Thompson scattering spectra, or 
by the measurement of Doppler widths of lines with negligible Stark widths. If then the 
discrepancy with, e.g., CCC calculations should still persist, one would have to accept the 
possibility of some broadening mechanism not yet considered. Finding any substantial error 
in the quantum calculations for conditions well inside the regime of validity of the impact 
approximation is not a realistic option, since they agree very well with, e.g., benchmark 
experiments of the most important inelastic cross section for the B III 2s-2p lines, namely 
the 2s-2p excitation cross section. 60 Although the elastic scattering contribution remains 
less accurate, even a 50% error in this contribution would increase the calculated width only 
by 15%. (For the Ne VII 3s-3p lines, such correction would require an increase of the elastic 
scattering contribution by a factor 2.5.) 

Although numerical improvements in the calculation of elastic scattering amplitudes 
would be very desirable, the pattern appears to be that the corresponding contribution to 
the width decreases with Z , while inelastic, dipole-allowed inelastic contributions dominate. 
They are quite accurately represented by the CBE approximation, which in turn explains 
the surprising accuracy of some semiempirical estimates. However, the perturbative or 
nonperturbative scmiclassical calculations cannot be trusted, if partial-wave contributions 
L £ 2-7T are important and/or if distances of closest approach approaching boundstate radii 
are responsible for much of the broadening. A practical criterion for the validity of these 
calculations is n 2 /Z ~ 2-7T. This appears to be consistent with our recent quantum-mechanical 
calculation for the NIV An = 0, n = 3 singlet line, which accounts for 70% of the measured 
width. For a more detailed analysis of the transition from the quantum-mechanical to the 
semiclassical regime, a good example is provided by an analogous discussion of electron 
scattering on bare nuclei. 61 
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